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ABSTRACT 

Deep X-band imaging of the most luminous z ~ 4 quasars currently offers the earliest pos- 
sible view of the mass-dominant stellar populations of the host galaxies which house the first 
super-massive black holes in the Universe. This is because, until the advent of the James Webb 
Space Telescope, it is not possible to obtain the necessary deep, sub-arcsec resolution imaging 
at rest-frame wavelengths Xrest > 4000A at any higher redshift. We here present and analyse 
the deepest, high-quality Ks-hand images ever obtained of luminous quasars at z ~ 4, in 
an attempt to determine the basic properties of their host galaxies less than 1 Gyr after the 
first recorded appearance of black holes with Mbh > 10^ Mq. To maximise the robustness 
of our results we have carefully selected two Sloan Digital Sky Survey quasars at z ~ 4. 
With absolute magnitudes A/,; < —28, these quasars are representative of the most luminous 
quasars known at this epoch but they also, crucially, lie within 40 arcsec of comparably-bright 
foreground stars (required for accurate point-spread-function definition), and have redshifts 
which ensure line-free Ks-hand imaging. The data were obtained in excellent seeing con- 
ditions (< 0.4-arcsec) at the European Southern Observatory on the Very Large Telescope 
with integration times of ~ 5.5 hours per source. Via carefully-controlled separation of host- 
galaxy and nuclear light, we estimate the luminosities and stellar masses of the host galaxies, 
and set constraints on their half-light radii. The apparent Ks-hand magnitudes of the quasar 
host galaxies are consistent with those of luminous radio galaxies at comparable redshifts, 
suggesting that these quasar hosts are also among the most massive galaxies in existence at 
this epoch. However, the quasar hosts are a factor ~ 5 smaller ((^1/2) = 1.8 kpc) than the 
host galaxies of luminous low-redshift quasars. We estimate the stellar masses of the z ~ 4 
host galaxies to lie in the range 2 — 10 x 10^^ M©, and use the CIV emission line in the 
Sloan optical spectra to estimate the masses of their central supermassive black holes. The 
results imply a black-hole:host-galaxy mass ratio Mbh : Mgai — 0.01 — 0.05. This is an 
order of magnitude higher than typically seen in the low-redshift Universe, and is consistent 
with existing evidence for a systematic growth in this mass ratio with increasing redshift (i.e. 
Mbh ■ Mgai ex (1 + z)^''^^"), at least for objects selected as powerful active galactic nuclei. 

Key words: galaxies: active - galaxies: fundamental parameters - galaxies: high-redshift - 
galaxies: photometry - infrared: galaxies - quasars: general. 



1 INTRODUCTION 

At all but the most modest redshifts, the measurement of the ve- 
locity width of broad permitted lines (e.g. H/3, Mgll, CIV) from 
unobscured Type-1 active galactic nuclei (AGN) offers the only di- 
rect way to estimate the masses of super-massive black holes (Wan- 
del, Peterson & Malkan 1999; Kaspi et al. 2000; McLure & Dunlop 
2001 ; Ho 2002; McLure & Jarvis 2002; Willott et al. 2003; McLure 
& Dunlop 2004; Vestergaard & Peterson 2006; Peng et al. 2006a; 
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Vestergaard et al. 2008; Willott et al. 2010). For this reason, de- 
termining the properties of AGN host galaxies remains of crucial 
importance for exploring the relationship between black-hole and 
galaxy evolution over cosmic history, and in particular for testing 
alternative theories for the origin of the now well-established pro- 
portionality between black-hole and galaxy bulge mass found in the 
local/present-day Universe (Kormendy & Richstone 1995; Magor- 
rian et al. 1998; Silk & Rees 1998; Gebhardt et al. 2000; Merritt 
& Ferrarese 2001; McLure & Dunlop 2002; Tremaine et al. 2002; 
Bettoni et al. 2003; Marconi & Hunt 2003, Merloni et al. 2010; 
Jahnke&Maccio2011). 

Quasars, as the most extreme of AGN, offer potentially the 
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greatest challenge to theory, especially at high redshift where, for 
example, black holes with masses Mbh > 1 x 10^ M© already 
appear to have been in place less than 1 billion years after the Big 
Bang (Fan et al. 2001, 2003; Willott et al. 2007; Mortlock et al. 
2011). However, extracting reliable information on the host galax- 
ies of these overwhelmingly bright AGN presents a formidable 
technical challenge. 

At low redshift {z < 0.4) this challenge was effectively 
met when the first refurbishment of the Hubble Space Tele- 
scope (HST) enabled optical imaging of the necessary depth and 
exquisite+stable angular resolution to allow the scalelengths and 
morphologies of low-redshift quasar hosts to be reliably determined 
for the first time (e.g. Disney et al. 1995; McLure et al. 1999; Dun- 
lop et al. 2003; Floyd et al. 2004). Indeed, the results of these deep 
HST optical imaging programs of quasar hosts, when combined 
with "virial" measurements of the black-hole masses in the same 
objects, were instrumental in demonstrating that these most active 
objects displayed the same black-hole:host-galaxy mass propor- 
tionality as found in quiescent galaxies (McLure & Dunlop 2002). 

However, extending the effective study of quasar host galaxies 
out to higher redshifts, and in particular into the "quasar epoch" at 
z > 2, has proved to be extremely difficult. There are a number of 
reasons for this. First, there is the obvious cosmological dimming 
of the host-galaxy light. Second there is the necessity of observ- 
ing in the near-infrared, in order to continue to sample the rest- 
frame optical light of the host galaxy (both for comparison with 
low-redshift studies, and to avoid the host-galaxy light being com- 
pletely swamped by the UV-bright quasar). Some success in detect- 
ing the hosts of quasars out to 2: ~ 2 was achieved with the NIC- 
MOS infrared camera on HST (Kukula et al. 2001, Ridgway et al. 
2001; Peng et al. 2006b) but, even with this high-resolution near- 
infrared imaging, obtaining reliable measurements of host-galaxy 
luminosities and scalelengths proved extremely difficult. We now 
know that this was almost certainly in part due to the fact that, 
in general, massive galaxies at high-redshift have since been dis- 
covered to be generally much more compact than their low-redshift 
counterparts (Daddi et al. 2005; Trujillo et al. 2006,2007; Longhetti 
et al. 2007; Zirm et al. 2007; Cimatti et al. 2008; van Dokkum et 
al. 2008; Buitrago et al. 2008; Szomoru et al. 2010; Targett et al. 
2011). If this is also true for quasar hosts, then the separation of 
galaxy and nuclear light in quasars ai z > 2 will inevitably be even 
more problematic than was perhaps first anticipated. 

Nonetheless, these difficulties have not deterred several 
groups from attempting to repeat this experiment at redshifts as 
high as z ~ 4 — 5. Such studies have been encouraged (at least in 
part) by the advent of active/adaptive optics on ground-based 8-m 
class telescopes, and by the fact that z ~ 4 is the natural limit for 
ground-based studies of quasar host galaxies if one hopes to de- 
tect host-galaxy light longward of the 4000A/Balmer break in the 
K-bsLnd. 

However, perhaps unsurprisingly, existing studies of the host 
galaxies of quasars at z ~ 4 have met with very limited suc- 
cess. For example, using HST NICMOS H-hand data, Peng et al. 
(2006b) attempted to measure the luminosities and sizes of the 
hosts of two gravitationally-lensed quasars at z — 4.1 and z = 4.5 
via two-dimensional modelling. Although reliable host-galaxy lu- 
minosities were claimed, the HST images sampled light shortward 
of Xrest ~ 4000 A and, in addition, the complexity of the lensing 
model led the authors to conclude that the extracted morphological 
parameters could not be trusted. In an alternative ground-based ap- 
proach, Hutchings (2003, 2005) observed seven quasars at z ~ 5 
with the Gemini 8-m telescope in the J, H, and ii'-bands. Using 



point-spread function (PSF) subtraction, host-galaxy luminosities 
were estimated from the PSF-subtracted images, although no ev- 
idence of smooth centrally-concentrated host galaxies was found 
and, again, this time due to the decision to observe quasars at z ~ 5, 
this imaging sampled rest-frame wavelengths Xrest < 4000 A, lim- 
iting the usefulness of the extracted luminosities. Most recently, 
McLeod & Bechtold (2009) observed 34 z ~ 4 quasars in the K- 
band with the Magellan I and Gemini North telescopes. However, 
despite the use of well-controlled PSF subtraction, only four host 
galaxies were detected in their deepest and sharpest data, at an av- 
erage quasar redshift (z) — 3.8. Host-galaxy luminosity and size 
estimates were attempted from the PSF-subtracted images of these 
4 quasars, although the authors noted that the surface-brightness 
limit of these data was really of inadequate depth for the accu- 
rate recovery of host-galaxy sizes. The deepest imaging obtained 
by these authors (~ 3.5 hours per source on 8-m class telescopes) 
would thus seem to represent a minimum requirement for the dis- 
covery and study of quasar host galaxies at z ~ 4. 

In this paper we present and analyse even deeper (~ 5.5 
hours per source) European Southern Observatory (ESO) 8.2-m 
Very Large Telescope (VLT) Ks-band imaging of two z ~ 4 Sloan 
Digital Sky Survey (SDSS) quasars. These two quasars were care- 
fully chosen to have a redshift which ensures line-free imaging of 
the host-galaxy stellar population at Xrest > 4000 A, and to have a 
star of comparable brightness within an angular radius of 40 arcsec, 
in order to provide a robust, high signaknoise ratio, on-detector rep- 
resentation of the VLT Ks -band PSF over the complete duration of 
the host-galaxy imaging. This care in quasar target selection, com- 
bined with deliberate (flexibly scheduled) use of only the very best 
seeing conditions (and the deliberate avoidance of adaptive optics), 
has enabled us to achieve robust detections of the host galaxies of 
both quasars, despite the fact that, with Mi < —28, these objects 
are representative of the most luminous/massive black holes in ex- 
istence at this (or indeed any) epoch. 

This paper is structured as follows. In Section 2 we describe 
our observing strategy, our target selection, and summarise the new 
near-infrared data obtained with the VLT. Next, in Section 3 we 
explain the PSF subtraction technique used to separate host and 
nuclear light. Then, in Section 4 we describe the determination 
of host-galaxy properties and the estimation of the central black- 
hole masses. Finally, in Section 5 we attempt to place our results in 
the context of studies of other galaxy populations at both high and 
low redshift, and briefly explore the consequences of this study for 
the inferred cosmic evolution of the black-hole:galaxy mass ratio. 
Throughout we adopt a cosmology with Ho — 70kms^^ Mpc^^, 
fim = 0.3, and SIa ~ 0.7. Unless otherwise noted, we report all 
magnitudes in the Vega system to ease comparison with previous 
studies, but our ii's-band photometry can be simply converted to 
the AB magnitude system via Ks,ab ~ Ks,vega + 1.85. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Instrument 

The deep Ks -band imaging presented here was performed with 
the 8.2-m VLT (located on Cerro Paranal in the Atacama Desert, 
northern Chile), using the Infrared Spectrometer And Array Cam- 
era (ISAAC), a 1024 x 1024 Hawaii Rockwell array with a pixel 
scale of 0.15 arcsec. This pixel scale provided the necessary good 
sampling of even the very best seeing disc, while still providing 
sufficient field-of-view to include the required nearby field stars for 
accurate on-image PSF characterisation. 
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2.2 Observational strategy 

The reliable detection of host-galaxy light around luminous quasars 
is difficult at all redshifts, but especially at high redshift, due to 
the small angular size of any reasonable anticipated host galaxy, 
the faintness of any potential host galaxy relative to the scattered 
nuclear light in the wings of the PSF, and the fact that the rest- 
frame optical light has been redshifted to the near-infrared. 

Predicting surface brightnesses (and hence exposure times) for 
high-2 quasar host galaxies is complicated by unknown amounts of 
stellar and morphological evolution to offset against the cosmologi- 
cal dimming and anticipated fc-corrections. In an attempt to best de- 
fine our observational strategy we generated synthetic VLT/ISAAC 
exposures of a z = 4 quasar with a nuclear/host ratio of 15, 0.6- 
arcsec seeing, and background-limited noise. The synthetic quasar 
had a total i^T^-band apparent magnitude of Ks = 16.5 (consis- 
tent with the anticipated A's-band magnitudes of our target SDSS 
quasars), comprising a nuclear point source with Ks = 16.8 and 
a de Vaucouleurs host galaxy with Ks ~ 19.7 (consistent with 
existing studies of the K — z relation for massive galaxies). For 
the simulated host galaxy we assumed a half-light radius of 5 kpc. 
Separate PSFs were used to construct the synthetic quasar and to 
carry out nuclear/host separation, with both PSFs obtained from 
the publicly-available ISAAC image of the GOODS-South field, 
on which they were ~ 30-arcsec apart (thus accurately reflecting 
our observing strategy of using a nearby star to characterise the on- 
image PSF). In order to separate host and nuclear components and 
reliably determine total luminosities for such galaxies, we found 
that we needed to follow the surface brightness profile out to ^^^ 3 
galaxy half-light radii. This required integrating down to a 3a sur- 
face brightness level of ^k — 25 mag arcsec^^. The VLT ISAAC 
Exposure Time Calculator indicated that we could achieve this level 
in ~ 5 hours on source, and so this was set as the minimum on- 
source exposure time for our observations. 



2.3 Target selection 

For the reasons explained above, we wished to obtain emission- 
line free images of the most luminous quasars at 2 ~ 4. Con- 
tamination of the observed if-band light by [OIII] and Hfi line 
emission is avoided by insisting on z > 3.95, and so we aimed to 
select quasar targets at only slightly higher redshift to ensure that 
the Ks -hand filter was still essentially filled by the continuum light 
starting at Xrest — 4000A. We confined our attention to the most 
luminous quasars in the SDSS quasar catalogue at these redshifts, 
corresponding to absolute magnitudes Mi < —28. These objects 
are rare; the final SDSS quasar catalogue contains only ~ 200 such 
quasars in the redshift range 4 < 2: < 5 over the full ~ 9380 deg^ 
of the SDSS quasar survey (Schneider et al. 2010). They thus have 
a comoving number density of ~ 2 per comoving Gpc"* at 2 ~ 4, 
which corresponds to ~ 1 object in the cosmological volume sim- 
ulated in the "Massive Black" simulation recently discussed by Di 
Matteo et al. (2011) (see Section 5). 

Finally, we further limited our observational options to 
quasars from within this sample which happened to have at least 
one bright star within a radius of 40 arcsec, in order to provide a 
high signahnoise ratio representation of the real form of the point- 
spread-function (PSF) on the final science image. It was deemed 
essential that the nearby stars were at least as bright as the quasar 
at Ks-band, but not so bright as to risk saturation or non-linearity 
in the background limited imaging. This final rather stringent re- 
quirement reduced our potential target list to less than 10 quasars. 



Table 1. Basic details of the VLT ISAAC A'g-band observations of the two 
z ~ 4 SDSS luminous quasar targets. Column 1 lists the target names. 
Column 2 gives the total exposure times in hours. Column 3 provides the 
measured 4-arcsec diameter aperture A'5-band magnitude (Vega) for each 
quasar Finally, Column 4 lists the spectroscopic redshifts of the quasars as 
given in the final SDSS quasar catalogue (Schneider et al. 2010). 



Source 



Exp. time / hr Ks (Vega) 



SDSS J131052.51-005533.2 
SDSS J144617.35-010131.1 



5.9 
5.1 



16.50±0.04 
16.70±0.05 



4.1567 
4.1608 



2.4 Data 

Our insistence on long (> 5 hr) integrations in the very best seeing 
conditions meant that we only obtained the required deep imaging 
for two of the potential target quasars, but all data were taken in 
photometric conditions and consistent high-quality (< 0. 4-arcsec) 
atmospheric seeing. Our observations of these two 2 ~ 4 quasars 
are summarised in Table 1. 



2.5 Data reduction 

Data reduction was performed using standard IRAF packages. Dark 
frames of equal integration time to the science data were taken on 
all nights. Each dark-subtracted image was divided by a normalised 
flat-field, derived though the scaled sigma-clipped median combi- 
nation of neighbouring science frames observed within 20 minutes 
of each image. Registering of the offset frames was performed us- 
ing the brightest stars in the images. A map of the bad pixels in 
the array was obtained and used to exclude these pixels during im- 
age combination. Median filtering was used to reject cosmic rays 
and produce an initial reduced image which was then processed 
to create a mask of all source flux detected in the mosaiced data. 
The science frames were then reprocessed to create an improved 
flat-field, and then the data were re-reduced to produce the final 
science sub-images. As accurate matching of both the quasar and 
PSF image quality is of critical importance in this work, separate 
mosaics were created for the quasars and for the PSF stars using 
offset shifts derived from each object individually. Final image reg- 
istration was performed using a SplineS interpolation in the IRAF 
routine IMSHIFT 



3 PSF SUBTRACTION 

Because of the expected overwhelming dominance of the nu- 
clear light, separation of nuclear and host-galaxy light was 
first performed using scaled PSF subtraction. We subtracted an 
appropriately-scaled version of the on-image stellar PSF from the 
quasar image to produce the minimum residual which still displays 
a monotonically declining surface-brightness profile (see Fig. 1), 
thus providing a robust estimate of the minimum level of host- 
galaxy flux. This method clearly represents an over-subtraction 
of the nuclear PSF, but at least ensures exclusion of all nuclear 
light from the subtracted residual. The likely extent of any over- 
subtraction of host-galaxy light is quantified via the host-galaxy 
modelling presented in Section 4. 

The Ks-band images of the two quasars, PSFs, and mono- 
tonic PSF-subtracted residuals are shown in Fig. 2. Several fac- 
tors complicate determination of the exact level of PSF to subtract. 
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Figure 1. Raw quasar (top) and PSF-subtracted monotonically-declining host-galaxy (bottom) .Kg-band surface brightness profiles for SDSS 
J131052.51-005533.2 (left) and SDSS J144617.35-010131.1 (right). 



where small differences between the shape of quasar and stellar 
PSFs can create an artificial residual. In order to check that the 
quasar-PSF subtracted residuals represent host-galaxy flux and not 
such anomalies, PSF-PSF subtractions were performed using alter- 
nate stars on the science images, with the scaling of the subtracted 
PSF again adjusted to the point where the residual (just) displayed 
a monotonically decreasing radial surface brightness profile. Re- 
assuringly, the resulting residuals from these tests (an example of 
which is shown in Fig. 2) were always much smaller than (and qual- 
itatively very different from) the putative quasar host galaxies. The 
results of these tests were also used to determine the typical pho- 
tometric error introduced by the small but inevitable mis-matches 
between the PSFs at different points in the ISAAC image. Finally, 
there is also inevitably some uncertainty when determining the pre- 
cise level of PSF normalisation required for subtraction to yield the 
minimum monotonically declining residual, with a narrow range of 
scaled PSFs producing acceptable monotonic residuals. 



4 RESULTS 

4.1 Host-galaxy luminosities 

In Table 2 we give the measured /i'g-band aperture magnitudes of 
the residual host galaxies, along with our estimates of the pho- 
tometric error contributions from basic uncertainty in the back- 
ground, PSF mis-match, and uncertainty in the PSF scaling re- 
quired to produce the monotonically declining residual image. 

In an attempt to quantify the extent to which the host galaxy 
luminosities will have been (inevitably) under-estimated by our 
adopted process of scaled PSF subtraction, we undertook two- 
dimensional modelling of the host galaxy residuals shown in Fig. 2 



using GALFIT (Peng et al. 2002, 2010). This modelling was also 
aimed at extracting information on host-galaxy size and morphol- 
ogy, as discussed further below. 

To minimise the impact on this modelling of the known over- 
subtraction of the PSF (which of course produces somewhat "flat- 
topped" residual images), the central 9 pixels of each host galaxy 
residual were excluded from the model fitting. The results of this 
model fitting are illustrated in Fig. 3, tabulated in Table 3, and dis- 
cussed further below 

From the point of view of host galaxy luminosities, the key 
point is that, as summarised in Table 3, the modelled restoration 
of additional light in the central regions of the host galaxies does, 
as expected, make the inferred host galaxies brighter, but only by 
~ 0.2 — 0.3 mag. within our adopted 4-arcsec diameter aperture. 



4.2 Host-galaxy sizes 

In modelling the quasar host-galaxies with GALFIT we allowed 
axial ratio, luminosity, scalelength (ri/2) and Sersic index (n) to 
vary. The resulting best-fitting parameter values are summarised 
in Table 3, and the resulting model images shown in Fig. 3. The 
derived best-fitting values of half-light radius and Sersic index (n) 
appear perfectly sensible, giving some additional reassurance that 
we are modelling genuine host-galaxy light. However, we found 
that while the half-light radius of the host-galaxies appears to be 
relatively robust, the derived value of Sersic index is not. 

This point is illustrated in Fig. 4 which shows how reduced 
xS changes as half-light radius is varied, for different alternative 
fixed values of Sersic index n. It can be seen that minimum xl is 
achieved at essentially the same length-scale for each of the three 
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Figure 2. The Ks-band quasar images, and the results of scaled quasar-PSF subtraction to produce residuals with monotonically declining radial luminosity 
profiles. The top row shows the image of the quasar J13 1052. 5 1—005533.2 (left), the stellar PSF extracted from within the same image (centre) and the 
PSF-subtracted monotonically declining host-galaxy residual (right). The second row shows the same information for the quasar SDSS J144617. 35— 010131.1 
with its associated on-image stellar PSF. The bottom row shows the effect of performing a control experiment, applying the same analysis to two stellar PSFs 
drawn from the same science image (left and centre) with the resulting residual from scaled PSF-PSF subtraction, again designed to produce a monotonically- 
declining residual (right). All panels are 12 X 12 arcsec. 

Table 2. Kg-hand photometry of the quasar host galaxies as derived from the scaled PSF-subtracted residuals. Column 1 gives the source name. Column 2 
gives the measured 4-arcsec diameter aperture magnitude with associated total error Columns 3, 4, and 5 list the contributions to this total error provided by 
uncertainty in sky background, from PSF mismatch, and from uncertainty in the precise level of PSF subtraction required to produce the desired minimum 
residual image with a monotonically declining radial luminosity profile. 



Source 



Ks I mag 



Sk back 5 k PSF <5r- sub 



SDSS J131052.51-005533.2 
SDSSJ144617.35-010131.1 



19.42±0.17 
20.51 ±0.29 



0.07 
0.12 



0.15 
0.26 



0.03 
0.07 
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Figure 3. Two-dimensional modelling of the host galaxies of SDSS J131()52.51-0()5533.2 (upper row) and SDSS J144617.35-010131.1 (lower row). The 
left-hand panel shows the host-galaxy light which remains after subtraction of the scaled PSF to produce the monotonically declining residual. The middle 
panel shows the best-fitting two-dimensional model after fitting to the residual host galaxy (excluding the central pixels). The right-hand panel shows the final 
residual after subtraction of the best fitting two-dimensional model from the observed host galaxy. As in Fig. 2, all panels are 12 X 12 arcsec. 



alternative values of n adopted (n = 2.5, 1, 4), but that the best 
quality of fit achieved is essentially independent of n. 

We conclude that we can say nothing about whether the host 
galaxies of these luminous 2 ~ 4 quasars are dominated by discs 
or spheroids, but are reassured that our data are of sufficient quality 
to constrain their scalelengths and luminosities. 



4.3 Host-galaxy stellar masses 

We have derived estimates of the stellar masses of the two quasar 
host galaxies using the evolutionary synthesis models of Bruzual 
& Chariot (2003), assuming the Initial Mass Function (IMF) of 
Salpeter (1955). The extent to which an accurate stellar mass can 
be derived is limited by the fact that, while we do at least know 
the precise spectroscopic redshift of both quasars, we possess no 
meaningful information on the colour of either host galaxy. We 
thus made two simple but extreme estimates of the stellar mass of 
each host galaxy (effectively bracketing the plausible mass range) 
by assuming first that the entire galaxy was formed as recently as 
Zf — 5, and then alternatively assuming that the galaxy was formed 
in a single burst at much higher redshift, Zf = 10. While neither 
assumption is completely realistic, this does at least yield fairly 
robust upper and lower limits on the possible stellar masses of the 
two quasar host galaxies. The resulting inferred minimum and max- 
imum stellar masses for each quasar host galaxy are given in Table 
4. We note that adoption of the IMF of Chabrier (2003) would yield 
stellar masses a factor ~ 1.8 lower. 



4.4 Black-hole masses 

The masses of the black holes which power these two z ~ 4 
quasars were estimated using the so-called virial black hole mass 
estimator for broad-line AGN (e.g. Wandel, Peterson & Malkan 
1999; Kaspi et al. 2000). Under the assumption that the broad-line 
emitting gas is in virial motion within the central black hole's grav- 
itational potential, the central mass can be estimated via Mbh ~ 
G^^RbirV^, where Rbir is the radius of the broad-line region 
(BLR) and V is the orbital velocity of the line-emitting gas. In 
practice, Rtir is estimated via the correlation between optical lu- 
minosity and Rbir discovered from reverberation mapping of low- 
redshift AGN (Kaspi et al. 2000), and the gas orbital velocity is 
taken to be the full width at half-maximum (FWHM) of one or 
more of the H/n, Mgll or CIV emission lines, depending on which 
lines are observable given the source redshift. 

At 2: ~ 4, the CIV line is the only one available within 
the SDSS optical spectra, and can be used to estimate Mbh using 
the virial black hole mass estimator from Vestergaard & Peterson 
(2006): 



Mbh 

'Me 



ALi 



1044 ergs- 



/' FWHM{CIV) \' 
V lOOOkms^i / 



-F6.66 



The results of this calculation for the two quasars considered 
here are given in Table 4. 
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Table 3. The results of two-dimensional modelling of the residual host-galaxy images. Column 1 gives the source name. Column 2 lists the model semi-major 
axis scalelength (half-light radius) in kiloparsec. Column 3 lists the value of the Sersic index (n). Column 4 gives the total model Kg-bmid magnitude. Column 
5 gives the model i^g-band magnitude within a 4-arcsec diameter aperture (for ease of comparison with the raw residual photometry values given in Table 2). 
Column 6 gives xi for the best-fitting model. 



Source 


ri/2/kpc 


Sersic n 


Ks model tot 


Ks model 4" 


xi 


SDSSJ131052.51-()05533.2 
SDSSJ144617.35-()1()131.1 


0.75 
2.84 


2.4 
2.8 


19.27 
20.24 


19.30 
20.29 


1.24 
1.29 




1.0 1.5 2.0 

half-light rodlus (kpc 



2 3 4 
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Figure 4. xl plotted against model half-light radius for SDSS J131052.51-005533.2 (left) and SDSS J144617.35-010131.1 (right) for different adopted 
morphologies as parameterised by Sersic index n. Results are shown for fixed values of Sersic index n = 1 (filled stars) and n = 4 (filled circles), as well as 
for the average best-fitting value of n = 2.5 (solid line). 



5 DISCUSSION 

In Fig. 5 (adapted from Miley & De Breuck 2008 following Rocca- 
Volmerange et al. 2004) we attempt to illustrate how the A's-band 
brightness of our two z ~ 4 quasar host galaxies compares both 
with other estimates of quasar host-galaxy brightness at compara- 
ble redshift, and with the brightness of other well-studied popula- 
tions of high-redshift galaxies. Here it can be seen that our quasar 
hosts lie on the now well-established K ~ z relation for powerful 
radio galaxies (e.g. Lilly & Longair 1984; Bales et al. 1997; van 
Breugel et al. 1998; De Breuck et al. 2002). As can be seen from 
Fig. 5, and as inferred from a series of deep near-infrared imaging 
studies of powerful radio galaxies (e.g. Dunlop & Peacock 1993; 
Best et al. 1998; Willott et al. 2003; Targett et al. 201 1), the radio- 
galaxy K — z relation essentially defines the high-mass envelope of 
the evolving galaxy population. Regardless of their precise stellar 
masses, it is thus clear that, as perhaps expected, the host galax- 
ies of the most luminous SDSS quasars at z ~ 4 are amongst the 
most massive known galaxies at this epoch. In addition, while we 
have argued that our data are of higher quality, and our results more 
robust than those derived from previous studies of quasar hosts at 
z ~ 4, our host galaxy luminosities are basically consistent with 
those derived by McLeod & Bechtold (2009), Hutchings (2003, 
2005) and Peng et al. (2006b). As judged from the galaxy evolu- 
tion models of Rocca-Volmerage et al. (2004) over-plotted on Fig. 
5, it would seem that all of the z ~ 4 quasar host galaxies uncov- 
ered to date have stellar masses in the range 1 — 10 x 10^^ M©, 
which is entirely consistent with our own derived bounds on the 
possible stellar masses of our two quasar host galaxies as given in 
Table 4. 

Our study is, however, the first to yield apparently robust val- 



ues for the half-light radii of quasar host galaxies at such high red- 
shifts. Despite their large stellar masses, the derived half-light radii 
of the 2: ~ 4 quasar hosts (ri/2 — 1.8 kpc) are clearly much smaller 
than those of either low-redshift quasar hosts (~ 10 kpc; Dunlop 
et al. 2003) or indeed the vast majority of present-day galaxies 
of comparable mass (e.g. Hyde & Bernardi 2009). However, our 
small derived values for ri/2 seem broadly as expected given the 
rapidly growing observational evidence for the increasing compact- 
ness of massive galaxies with increasing redshift (Daddi et al. 2005 
Trujillo et al. 2006,2007; Longhetti et al. 2007; Zirm et al. 2007 
Cimatti et al. 2008; Buitrago et al. 2008; van Dokkum et al. 2008 
Targett et al. 201 1) and the expectations from recent simulations of 
galaxy growth (e.g. Oser et al. 201 1). We conclude that, just as at 
low redshift, the properties of quasar host galaxies at z ~ 4 seem 
broadly as expected for "normal" galaxies of comparable mass at 
that particular cosmological epoch. 

The derived black-hole masses given in Table 4 are clearly 
very large, approaching M^h — W^^ M© for both quasars, and 
so it is worth pausing to consider whether such values are plausi- 
ble/expected. Our conclusion is that these masses, while high, are 
certainly not unreasonable given the relative rarity of these ultra- 
bright quasars, our understanding of quasar demographics, and the 
results of the latest high-redshift simulations. Specifically, our tar- 
get objects were selected from a parent sample which contains 
only 1 such luminous quasar in the redshift range 4 < z < 5 
per 50 deg^ on the sky, corresponding to a comoving source den- 
sity of only ~ 2Gpc^^. These are thus extreme objects, ar- 
guably the natural descendents of the most extreme quasars found 
at z ~ 6.5 — 7, which already appear to have black-hole masses 
Mth ~ 2 X 10^ M0 (Willott, McLure & Jarvis 2003; Mortlock 
et al. 2011). Only 60 million years of further Eddington-limited 
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Table 4. Estimated galaxy stellar masses and quasar black-hole masses. Column 1 gives the source name. Column 2 gives the host-galaxy stellar mass (in units 
of 10^ Mq) on the assumption that the entire galaxy was formed in a single starburst at zj = 5, while column 2 gives the largest value of stellar mass inferred 



from pushing the formation epoch back to Zf 
Section 4.4. 



10. Column 4 gives the value of black-hole mass estimated from the SDSS optical spectra as described in 



Source 
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SDSS 1131052.51-005533.2 
SDSS 1144617.35-010131.1 
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Figure 5. Composite K — z diagram showing 2 ~ 4 quasar hosts, radio galaxies, and optically-selected galaxies, adapted from Miley & De Breuck (2008) 
following Rocca-Volmerange et al. (2004). The z ~ 4 SDSS quasar host-galaxy magnitudes derived here are plotted (solid black circles) along with the other 
z ^ A quasar host-galaxy results (solid grey circles) from the literature (Peng et al. 2006b; Hutchings 2003, 2005; McLeod & Bechtold 2009). Also plotted are 
massive radio galaxies (squares) from De Breuck et al. (2002), and optically-selected field galaxies in the Hubble Deep Field North (crosses) and the Hawaii 
survey (diamonds). These data are over plotted with theoretical K — z tracks for elliptical galaxies formed with a range of assumed initial reservoir gas masses, 
and developed using Pegase galaxy evolution models as described in Rocca-Volmerange et al. (2004). 



growth would be required to boost these black-hole masses by a 
further factor of 4, to values consistent with those derived here 800 
million years later. This viewpoint receives some further theoretical 
support from recent simulations of black-hole growth in the young 
Universe, such as the "Massive Black" simulation of Di Matteo et 
al. (201 1). The luminous quasars targeted here are sufficiently rare 
that, on average, only one would be expected in the ~ 0.5 Gpc"' 
comoving volume followed in this simulation, and Di Matteo et al. 
report that the most massive black hole in their simulation has al- 
ready achieved a mass of Mth — 5 x 10^ Mq by 2 ~ 5. By this 
point feedback and shock heating of infalling gas are inferred to 
limit further growth, but clearly the concept that such rare black 
holes may achieve masses approaching Mi,h — 10^° Mq by 2 ~ 4 
is consistent with current observational and theoretical constraints. 

Finally, we can combine our estimates of black-hole mass and 
galaxy stellar mass to infer the black-hole:galaxy mass ratio for 
our z ~ 4 quasars. The averaged value of Mth '■ Mgai is shown 
in Fig. 6, where we plot two 2 ~ 4 data points to indicate the dif- 
ferent values which arise from adoption of the either the maximum 
or minimum possible galaxy stellar masses (as derived assuming a 



formation redshift of Zf = 10 or 2/ = 5 respectively). The light- 
grey shaded area illustrates the ±0.3-dex uncertainty on the local 
Mbh/Mbuige ratio, centred on Mbh/Mbnige = 0.002. The sohd 
line shows the best fit to the observed evolution in 3CRR galax- 
ies from McLure et al. (2006), while the dark-grey shaded area 
shows the la uncertainty of this fit. This evolving Mbh/Mbuige 
ratio result is consistent with results from Peng et al. (2006a), 
who find a similar evolution out to 2 ~ 2 in a sample of 11 
quasars. Also shown is the 2 = 6.41 quasar observed by Willott, 
McLure, & Jarvis (2003) who used iiT-band and K-banA spec- 
tra covering the Mgll broad emission line to derive a black hole 
mass of 3 X 10^ Mq. However, note that in this case the host- 
galaxy mass refers not to a stellar mass inferred from the sort of 
host-galaxy imaging attempted here, but is a dynamical mass in- 
ferred from Very Large Array measurements of molecular gas in 
the quasar host galaxy by Walter et al. (2004). Even given the 
uncertainties, our results imply a black-hole:host-galaxy mass ra- 
tio (Mbh '■ Mgal — 0.01 — 0.05) which appears to be a fac- 
tor of ~ 10 higher than typically seen in the low-redshift Uni- 
verse. Fig. 6 demonstrates that this is consistent with the grow- 



© 0000 RAS, MNRAS 000, 000-000 



Host galaxies of luminous q 



4 - 



10 



Figure 6. Average Mf,^ : Mgai of the z ~ 4 SDSS quasars (filled circles) for both the maximum (lower point) and minimum (upper point) galaxy masses, 
derived for the host galaxies where formation redshift was constrained to be Zf = 10 or Zj- = 5 respectively. The solid line shows the best fit to the observed 
evolution in 3CRR galaxies (open circles) from McLure et al. (2006), while the dark-grey shaded area shows the Icr uncertainty of this fit. This trend is 
consistent with results from Peng et al. (2006b), who find a similar evolution out to z ~ 2 in a sample of 1 1 quasars (open triangle). The light-grey shaded area 
illustrates the ±0.3-dex uncertainty on the local Mf,f^/Mfj^ig^ ratio, centred on Mf^f^/Mfj^i^^ = 0.002. The data point at shown for the z = 6.41 studied by 
Willott, McLure, & Jarvis (2003) is based not on any estimate of the stellar mass of the host galaxy, but on a dynamical mass inferred from Very Large Array 
measurements of molecular gas in the quasar host galaxy by Walter et al. (2004). 



ing body of evidence for a systematic growth in this mass ratio 
with increasing redshift at least for objects selected as powerful 
active galactic nuclei. Indeed, an order of magnitude increase in 
Mbh '■ Mgai is exactly as anticipated given the redshift dependence 
of Mbh ■ Mgai oc (1 + z)^** inferred from recent studies of X-ray 
selected AGN (Bennert et al. 201 1). 



6 CONCLUSION 

We have obtained and analysed the deepest, high-quality ifg-band 
images ever obtained of luminous quasars at z ~ 4, in an at- 
tempt to determine the basic properties of their host galaxies less 
than 1 Gyr after the first recorded appearance of black holes with 
M^,h > 10^ Mq . To maximise the robustness of our results, we 
carefully selected two SDSS quasars at z ~ 4. With absolute mag- 
nitudes Mi < —28, these quasars are representative of the most 
luminous quasars known at this epoch but they also, crucially, lie 
within 40 arcsec of comparably-bright foreground stars (required 
for accurate point-spread-function definition), and have redshifts 
which ensure line-free Ks -band imaging. 

Our new data were obtained in excellent seeing conditions 
« 0.4-arcsec) with ISAAC on the ESO VLT, with integration 
times of ~ 5.5 hours per source. Via carefully-controlled separa- 
tion of host-galaxy and nuclear light, followed by two-dimensional 
galaxy model-fitting, we have estimated the luminosities and stellar 
masses of the host galaxies. In addition, although it did not prove 
possible to determine basic morphological type (via, e.g. Sersic in- 
dex n), we have extracted what appear to be robust estimates of the 
host-galaxy half-light radii. 

We find that the Ks-band magnitudes of the quasar host galax- 
ies are consistent with those of luminous radio galaxies at com- 



parable redshifts. This suggests that these quasar hosts are also 
among the most massive galaxies in existence at this epoch, and 
indeed inferred stellar masses of these host galaxies lie in the range 
2 — 10 X 10^^ Mq, depending on assumed formation redshift. We 
also find that these massive 2 ~ 4 quasar host galaxies are a factor 
~ 5 smaller ({ri/2) = 1.8 kpc) than the host galaxies of luminous 
low-redshift quasars, or indeed quiescent galaxies of comparable 
stellar mass. However, we argue that this is as expected given the 
growing evidence for the increasing compactness of massive galax- 
ies with increasing redshift, and conclude that, just as at low red- 
shift, the properties of quasar host galaxies at 2 ~ 4 are broadly as 
expected for "normal" galaxies of comparable mass at the relevant 
cosmological epoch. 



Finally we have used the CIV emission line in the SDSS opti- 
cal spectra of the quasars to estimate the masses of their central 
supermassive black holes. The results imply extreme black-hole 
masses approaching Alth — 10^'' Mq, but we argue that this is not 
unrealistic given the extreme nature and rarity of our target 2 ~ 4 
SDSS quasars (i.e. 1 — 2 per comoving Gpc'^), and current obser- 
vational and theoretical understanding of black-hole growth in the 
young Universe. 



Combining the black-hole and host-galaxy mass measure- 
ments we infer a black-hole:host-galaxy mass ratio Mbh ■ Mgai — 
0.01 — 0.05. This is an order of magnitude higher than typically 
seen in the low-redshift Universe, and is consistent with existing 
evidence for a systematic growth in this mass ratio with increas- 
ing redshift, at least for objects selected as powerful active galactic 
nuclei. 
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